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Q^- ABSTRACT 
D _ 

, Context. Many current and future surveys aim to detect the highest redshift (z > 7) sources through their Lyman-a (Lya) emission, 

using the narrow-band imaging method. However, to date the surveys have only yielded non-detections and upper limits as no survey 
has reached the necessary combination of depth and area to detect these very young star forming galaxies. 

Aims. We aim to calculate model luminosity functions and mock surveys of Lya emitters at 2 > 7 based on a variety of approaches 
, calibrated and tested on observational data at lower redshifts. 

^ |. Methods. We calculate model luminosity functions at different redshifts based on three different approaches: a semi-analytical model 

I ' based on CDM, a simple phenomenological model, and an extrapolation of observed Schechter functions at lower redshifts. The 

O ' results of the first two models are compared with observations made at redshifts z ~ 5.7 and z ~ 6.5, and they are then extrapolated 

to higher redshift. 

t/^ . Results. We present model luminosity functions for redshifts between z — 7 — 12.5 and give specific number predictions for future 

^ , ' planned or possible narrow-band surveys for Lya emitters. We also investigate what constraints future observations will be able to 

place on the Lya luminosity function at very high redshift. 

Conclusions. It should be possible to observe z — 7 ~ 10 Lya emitters with present or near-future instruments if enough observing 
^ ■ time is allocated. In particular, large area surveys such as ELVIS (Emission Line galaxies with VISTA Survey) will be useful in 

QQ ' collecting a large sample. However, to get a large enough sample to constrain well the z > 10 Lya luminosity function, instruments 

' further in the future, such as an ELT, will be necessary. 

' Key words, cosmology: theory - cosmology: early universe - galaxies: high redshift - surveys 

, 1 . Introduction but have only yielded upper limits. Future surveys planned for 



these redshifts include DaZle (Dark Ages Z Lyman-a Explorer, 
O ■ One of the most promising ways of detecting very high red- Norton et al. 2004) and ELVIS (Emission-Line galaxies with 
■ shift iz > 5), star-forming galaxies is via narrow-band imag- yiSTA Survey, Nilsson et al. 2006b). Observations of very high 
ing surveys targeting Lyman-a (Lya). In particular, redshifts jg^shift LEGOs have been proposed as an excellent probe of 
z ~ 5.7 and 6.5 have been extensively surveyed by several reionisation, through its effects on the Lya emission line pro- 
groups (e.g. Ajiki et al. 2003; Hu et al. 2004; Shimasaku et al. ^jg (g g Miralda-Escude 1998; Miralda-Escude & Rees 1998; 
2005; Ouchi et al. 2005, 2007; Malhotra et al. 2005; Tamguchi Caiman 2002; Gnedin & Prada 2004), the luminosity function 
et al. 2005; Tapken et al. 2006; Kashikawa et al. 2006). The cur- (g g Caiman & Gen 2005; Dijkstra, Wyithe & Haiman 2007) 
rent redshift record for a spectroscopically confirmed Lya emit- ^^^j clustering of sources (McQuinn et al. 2007). 
ter (LEGO - Lya Emitting Galaxy-building Object; see M0ller 

& Fynbo 2001) is 2 = 6.96 (lye et al. 2006) although Stark We here focus on Lya emission from star-forming galaxies, 
et al. (2007) have suggested the discovery of two LEGOs at where the Lya photons are emitted from gas which is photo- 
z = 8.99 and 9.32. The reason why narrow-band surveys are ionised by massive young stars. During recent years, theoretical 
restricted to a discrete number of narrow redshift windows is the work on Lya emitting galaxies has made significant progress, 
night sky OH emission lines. According to the OH line adas of There are three main aspects to these studies: i) predicting the 
Rousselot et al. (2000), at Lya redshifts ZLya > 7 (A > 9800 A) numbers of star-forming galaxies as a function of star formation 
there are only a few possible wavelengths where a narrow-band rate and redshift, ii) calculating the fraction of the Lya photons 
filter can fit in between the OH sky lines. These correspond to which escape from galaxies into the IGM and Hi) calculating the 
ZLya ~ 7.7, 8.2, 8.8, 9.4 and 10.1 - 10.5. Several future sur- factor by which the Lya flux is attenuated by scattering in the 
veys will target these windows in the sky aiming to detect very IGM on its way to the observer. Accurate treatments of ii) and Hi) 
high redshift galaxies. Three narrow-band surveys for Lya at are complicated because Lya photons are resonantly scattered 
redshift z ~ 8.8 have already been completed (Parkes, Collins by hydrogen atoms, with the consequences that absorption of 
& Joseph 1994; Wiflis & Courbin 2005; Cuby et al. 2007) Lya by dust in galaxies is hugely amplified, thereby reducing the 
escape fraction, and that even a small neutral fraction in the IGM 
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most theoretical papers have chosen to concentrate on only one 
aspect, adopting simplified treatments of the other two aspects. 
Haiman & Spaans(1999) made predictions of the number counts 
of Lya emitting galaxies by combining the Press-Schechter for- 
malism with a treatment of the inhomogeneous dust distribution 
inside galaxies. Barton et al. (2004) and Furlanetto et al. (2005) 
calculated the numbers of Lya emitters in cosmological hydro- 
dynamical simulations of galaxy formation, but did not directly 
calculate the radiative transfer of Lya photons. Radiative trans- 
fer calculations of the escape of Lya photons from galaxies in- 
clude those of Zheng & Miralda-Escude (2002), Ahn (2004) and 
Verhamme, Schaerer & Maselli (2006) for idealised geometries, 
and Tasitsiomi (2006) and Laursen & Sommer-Larsen (2007) 
for galaxies in cosmological hydrodynamical simulations. The 
transmission of Lya through the IGM has been investigated 
by Miralda-Escude (1998), Haiman (2002), Santos (2004) and 
Dijkstra, Lidz & Wyithe (2007), among others. Several authors 
(e.g. Haiman, Spaans & Quataert 2000; Fardal et al. 2001; 
Furlanetto et al. 2005) have studied the effect of cold accretion 
to describe the nature of so-called Lya blobs (Steidel et al. 2000; 
Matsuda et al. 2004; Nilsson et al. 2006a), see also sec.|6] 

Two models in particular, dissimilar in their physical as- 
sumptions, have been shown to be successful in reproducing 
the observed number counts and luminosity functions of Lya 
emitting galaxies at high redshifts: firstly, the phenomenological 
model of Thommes & Meisenheimer (2005) which assumes that 
Lya emitters are associated with the formation phase of galaxy 
spheroids, and secondly the semi-analytical model GALFORM 
(Cole et al. 2000, Baugh et al. 2005), which follows the growth 
of structures in a hierarchical, ACDM scenario. The GALFORM 
predictions for Lya emitters are described in detail Le Delliou 
et al. (2005, 2006) and Orsi et al. (in prep.), who show that the 
model is successful in reproducing both the luminosity functions 
of Lya emitting galaxies in the range 3 < z < 6 and also their 
clustering properties. 

In this paper we aim to provide model predictions to help 
guide the design of future planned or possible narrow-band sur- 
veys for very high redshift Lya emitters. We make predictions 
based on three approaches; the semi-analytical and phenomeno- 
logical models aheady mentioned, and an extrapolation from ob- 
servations at lower redshift. In section|2]we describe the different 
models used to make the predictions, and in section|3]we present 
the predicted number counts and comparisons with observed lu- 
minosity functions at lower redshifts. In section|4]we make num- 
ber predictions for some specific future surveys. A brief discus- 
sion regarding what can be learned from these future surveys is 
found in section |5] We give our conclusions in section |6] 

Throughout this paper, we assume a cosmology with Hq = 
70 km s^i Mpc"\ fim = 0.3 and ^Ia = 0.7, apart from the 
mock surveys discussed in section |5] which use GALFORM 
models matched to the cosmology of the Millenium Run 
(Springel et al. 2005), (which has iTg = 73 km s^^ Mpc"\ 

= 0.25 and = 0.75). 



2. Models 

We use three different approaches to predict the numbers of high 
redshift {z > 7) Lya emitters. The models are based on very dis- 
parate assumptions. The first model is the semi-analytical model 
GALFORM (Le Delliou et al. 2005, 2006), the second is the 
phenomenological model of Thommes & Meisenheimer (2005), 



and the third model is based on directly extrapolating from ob- 
servational data at lower redshifts. 

Both the semi-analytical and phenomenological models as- 
sume that the fraction of Lya photons escaping from galax- 
ies is constant, and that the IGM is transparent to Lya. The 
simple expectation is that before reionisation, the IGM will be 
highly opaque to Lya, and after reionisation it will be mostly 
transparent. However, various effects can modify this simple be- 
haviour; e.g. Santos (2004) finds that the transmitted fraction 
could be significant even before reionisation, while Dijkstra, 
Lidz & Wyithe (2007) argue that attenuation could be impor- 
tant even after most of the IGM has been reionised. The WMAP 
3-year data on the polarisation of the microwave background im- 
ply that reionisation occurred in the range z ~ 8 — 15 (Spergel et 
al. 2007), i.e. the IGM may be mostly transparent to Lya at the 
redshifts of most interest in this paper. In any case, what is im- 
portant for predicting fluxes of Lya emitters is the product of the 
escape fraction from galaxies with the attenuation by the IGM. 
The two effects are in this respect degenerate. 

2.1. Semi-analytical model 

The semi-analytical model GALFORM (Cole et al. 2000; Baugh 
et al. 2005), which is based on ACDM, has been shown to be suc- 
cessful in reproducing a range of galaxy properties at both high 
and low redshift, including Lya emitters in the range z = 3 — 6 
(Le DelHou et al. 2005; 2006). A full description of GALFORM 
is given in these earlier papers, so we only give a brief summary 
here. GALFORM calculates the build-up of dark halos by merg- 
ing, and the assembly of the baryonic mass of galaxies through 
both gas cooling in halos and galaxy mergers. It includes pre- 
scriptions for two modes of star formation - quiescent star for- 
mation in disks, and starbursts triggered by galaxy mergers - and 
also for feedback from supernovae and photo-ionisation. Finally, 
GALFORM includes chemical evolution of the gas and stars, 
and detailed stellar population synthesis to compute the stellar 
continuum luminosity from each galaxy consistent with its star 
formation history, IMF and metallicity (see Cole et al. 2000 for 
more details). The unextincted Lya luminosity of each model 
galaxy is then computed from the ionising luminosity of its stel- 
lar continuum, assuming that all ionising photons are absorbed 
by neutral gas in the galaxy, with case B recombination. 

The semi-analytical approach then allows us to obtain the 
properties of the Lya emission of galaxies and their abundances 
as a function of redshift, calculating the star formation histories 
for the entire galaxy population, following a hierarchical evolu- 
tion of the galaxy host haloes. In addition, when incorporated 
into an N-body simulation, we also obtain spatial clustering in- 
formation. This model has been incorporated into the largest N- 
body simulation to date, the Millennium Simulation (Springel et 
al. 2005), to predict clustering properties of Lya galaxies. These 
results will be presented in a forthcoming paper (Orsi et al., in 
prep.). 

The version of GALFORM which we use here is the one de- 
scribed in Baugh et al. (2005) and Le DeUiou et al. (2006), with 
the same values for parameters. The parameters in the model 
were chosen in order to match a range of properties of present- 
day galaxies, as well as the numbers of Lyman Break and sub- 
mm galaxies at z ~ 2 — 3. We assume a Kennicutt IMF for 
quiescent star formation, but a top-heavy IMF for starbursts, in 
order to reproduce the numbers of sub-mm galaxies. The only 
parameter which has been adjusted to match observations of Lya 
emitters is the Lya escape fraction, which is taken to have a con- 
stant value fesc = 0.02, regardless of galaxy dust properties. 
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Le Delliou et al. (2006) show that the simple choice of a con- 
stant escape fraction fesc = 0.02 predicts luminosity functions 
of Lya emitters in remarkably good agreement with observa- 
tional data at 3 < 2 < 6. Le Delliou et al. (2006) also compared 
the predicted Lya equivalent widths with observational data at 
3 < z < 5, including some model galaxies with rest-frame 
equivalent widths of several lOOA, and found broad consistency. 
For this reason, we use the same value f^sc ~ 0.02 for making 
most of our predictions at z > 7. However, since the value of the 
escape fraction at z > 7 is a priori uncertain in the models (e.g. 
it might increase with redshift if high redshift galaxies are less 
dusty) we also present some predictions for other values of fesc- 
Reionisation of the IGM affects predictions for the numbers 
of Lya emitters in deep surveys in two ways: i) feedback from 
photo-ionisation inhibits galaxy formation in low-mass halos 
and ii) reionisation changes the opacity of the IGM to Lya pho- 
tons travelling to us from a distant galaxy, as discussed above. 
GALFORM models the first effect in a simple way, approximat- 
ing reionisation as being instantaneous at redshift Zreion (see Le 
Delliou et al. 2006 for more details). We assume Zreion = 10, 
in Hne with the WMAP 3-year results (Spergel et al. 2007). 
As was shown in Le Delliou et al. (2006; see their Fig. 8), 
as far as the feedback effect is concerned, varying zreion over 
the range 7 < Zreion ^ 10 does not have much effect on 
the bright end of the Lya luminosity function most relevant to 
current and planned surveys. For example, varying Zreion be- 
tween 7 and 10 changes the predicted luminosity function at 
LLya > 10"^-^ erg s-1 by less than 10% for z - 7 - 10. 

2.2. Phenomenological model 

The phenomenological model of Thommes & Meisenheimer 
(2005; TM05 hereafter) assumes that the Lya emitters seen at 
high redshift are galaxy spheroids seen during their formation 
phase. We summarise the main features here, and refer the reader 
to TM05 for more details. The model is normalised to give the 
observed mass function of spheroids at z = 0, which is com- 
bined with a phenomenological function that gives the distri- 
bution of spheroid formation events in mass and redshift. Each 
galaxy is assumed to be visible as a Lya emitter during an initial 
starburst phase of fixed duration (and Gaussian in time), during 
which the peak SFR is proportional to the baryonic mass and in- 
versely proportional to the halo collapse time. The effects of the 
IMF and the escape fraction on the Lya luminosity of a galaxy 
are combined into a single constant factor (i.e. the escape frac- 
tion is effectively assumed to be constant). With these assump- 
tions, the luminosity function of Lya emitters can be computed 
as a function of redshift. The free parameters in the model were 
chosen by TM05 to match the observed number counts of Lya 
emitters at 3.5 < z < 5.7 (analogously to the choice of fesc in 
the GALFORM model). This model does not include any effects 
from reionisation. 

2.3. Observational extrapolation 

Our third approach is to assume that the Lya luminosity function 
is a Schechter function at all redshifts, following 

(t){L)dL = (j)*{L/L*)°'cxp{-L/L*)dL/L* (1) 

and to derive the Schechter parameters a, <j>* and L* at high 
redshifts by extrapolating from the observed values at lower red- 
shifts. For our extrapolation, we use fits to observations at red- 
shift z w 3 (van Breukelen et al. 2005; Gronwall et al. 2007; 



Table 1. Parameters of the fitted Schechter function in pre- 
viously published papers. References are 1) van Breukelen 
et al. (2005), 2) Gronwall et al. (2007), 3) Ouchi et al. 
(2007), 4) Dawson et al. (2007), 5) Malhotra & Rhoads 
(2004), 6) Shimasaku et al. (2006), and 7) Kashikawa et al. 
(2006). References 3 — 6 fit for three faint end slopes (a = 
— 1.0, —1.5 and —2.0), but here we only reproduce the results 
for fits with a — —1.5 as we fix the slope in our calculations. 
Malhotra & Rhoads (2004) do not give error bars on the fits. 
Dawson et al. (2007) fix the slope to a = —1.6. 



Ref 


Redshift 


a 


log (j!>*Mpc~^ 


log L*ergs/s 


1 


~ 3.2 


-1.6 




42.70^"o;i^ 


2 


3.1 


_1 40+0.45 
-'-•^-'-0.54 


-2.84 


42.46«:?^ 


3 


3.1 


-1.5 


-3.04«:i; 


42.76«:««6 


3 


3.7 


-1.5 




43.01«:«? 


4 


4.5 


-1.6 


q 77+0.05 
■J- ' ' -0.05 


43.04t°:ll 


5 


5.7 


-1.5 


-4.0 


43.0 


6 


5.7 


-1.5 


_n 44+0.20 


43.04t°-^i2 


3 


5.7 


-1.5 


q 1 1 +0.29 
■^•^^-0.31 


42.83t°:l« 


5 


6.5 


-1.5 


-3.3 


42.6 


7 


6.5 


-1.5 


26 


42.6ot;j:lf, 



Table 2. Extrapolated parameters of the observed Schechter 
function at higher redshifts. The faint end slope is fixed to 
a = -1.5. 



Redshift 


log <^*Mpc"-^ 


log L*crgs/s 


7.7 


-3.73 ±0.50 


42.88 ± 0.24 


8.2 


-3.80 ±0.50 


42.89 ± 0.24 


8.8 


-3.88 ±0.50 


42.91 ±0.24 


9.4 


-3.96 ±0.50 


42.92 ±0.24 


12.5 


-4.38 ±0.50 


42.99 ± 0.24 



Ouchi et al. 2007), z = 3.7 (Ouchi et al. 2007), z = 4.5 (Dawson 
et al. 2007), z w 5.7 (Malhotra & Rhoads 2004; Shimasaku et 
al. 2006; Ouchi et al. 2007) and z « 6.5 (Malhotra & Rhoads 
2004; Kashikawa et al. 2006), as found in Table [T] We make 
linear fits to log (f>* and log L* vs z, and extrapolate to higher 
redshift. For simplicity, we assume a fixed faint end slope of 
a = —1.5. We do not make any corrections for any possible ef- 
fects of reionisation or IGM opacity. The extrapolated values are 
given in Table ID 

3. Luminosity functions 

The possible Lya redshifts between z = 7 and z ~ 
10 where a narrow-band filter can be placed are z^ya = 
7.7, 8.2, 8.8, and 9.4. Redshifts beyond 10 are unreachable with 
ground-based instruments of the near-future. However, one pos- 
sibility for z > 10 surveys may be the James Webb Space 
Telescope (JWST, see section l43b and so we also make predic- 
tions for redshift z^ya = 12.5. 

First, we compare the Lya luminosity functions predicted 
by the semi-analytical (GALFORM) and phenomenological 
(TM05) models with current observational data at z ~ 6. This 
comparison is shown in Fig. [T] where we compare the models 
with the cumulative luminosity functions measured in several 
published surveys at z = 5.7 and z = 6.5. We can see that both 
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tweenredshiftS.S and 12.5, the Universe went through an impor- 
tant transition at this time as reionisation occurred (Spergel et al. 
2007). However, we do not know exactly how and when this hap- 
pened. Also, during this epoch the structure in the dark matter 
(and hence also in galaxies) was building up very rapidly. This 
underlines the interest of obtaining observational constraints at 
these redshifts. 

The hatched regions in Fig. |2] show the region of the lumi- 
nosity function diagram that has been observationally excluded 
at z = 8.8 by Willis & Courbin (2005) and Cuby et al. (2006). 
The former survey was deeper but in a smaller area, whereas the 
latter was more shallow over a larger area, hence the two-step 
appearance of the hatched area. From the plot, it is obvious that 
their non-detections are perfectly consistent with our theoretical 
models, although the GALFORM model with the non-standard 
escape fraction fesc ~ 0.2 is marginally excluded. 



Fig. 1. Plot of luminosity functions at redshifts z = 5.7 and 6.5. 
Red points and lines are at redshift z ~ 5.7, black points/lines at 
redshift z = 6.5. Points are observations by Ajiki et al. (2003; 
redshift 5.7, squares), Hu et al. (2004; redshift 5.7, pluses), 
Taniguchi et al. (2005; redshift 6.5, crosses), Shimasaku et al. 
(2006; redshift 5.7, diamonds), Kashikawa et al. (2006; red- 
shift 6.5, stars), Malhotra & Rhoads (2004; redshift 5.7 and 6.5, 
triangles) and Ouchi et al. (2007; redshift 5.7, filled squares). 
Solid lines show the GALFORM model (with escape fraction 
/gj^ = 0.02), dot-dashed Hnes the TM05 model. Note that the 
Taniguchi et al. (2005) and the Ouchi et al. (2007) samples are 
the spectroscopic samples only. 



4. Future surveys 

In this section, we discuss more specific predictions for several 
planned and possible future surveys. For all calculations, we as- 
sume a simple selection on the flux of the Lya emission line, 
with no additional selection on the equivalent width (i.e. we in- 
clude all galaxies with EWLya > 0). We also assume no ab- 
sorption by the neutral hydrogen in the IGM which would re- 
duce the measured fluxes and for GALFORM predictions we 
assume an escape fraction of fesc ~ 0.02. The predictions from 
the GALFORM and TM05 models for these future surveys as 
well as some published surveys are summarised in Table [3] 



models match the observational data reasonably well, once one 
takes account of the observational uncertainties. The error bars 
on the observational data points, omitted in the plot in order to 
not confuse the points, are large, at the bright end of the luminos- 
ity function due to small number statistics, and at the faint end 
due to incompleteness in the samples. The shallow slopes at the 
faint ends of the Taniguchi et al. (2005) and Ouchi et al. (2007) 
luminosity functions may be due to spectroscopic incomplete- 
ness. Both models fit the observations well. Hence we conclude 
that both of these models can be used to extrapolate to higher 
redshifts. 

We now have three methods of extrapolating to higher red- 
shifts, when the direct extrapolation of the Schechter function 
from lower redshifts is included. In Fig.|2]we plot the predicted 
luminosity functions at z = 7.7, 8.8 and 12.5 computed by these 
three methods. For other redshifts, the curves may be interpo- 
lated. For GALFORM, we show predictions for the standard 
value of the escape fraction fesc = 0.02 in the left panel, and 
for a larger value fesc = 0.2 in the right panel. This illustrates 
the sensitivity of the predictions to the assumed value of fesc 
at high redshift. The predictions from the other two models are 
plotted identically in both panels, since they do not explicitly 
include the escape fraction as a parameter GALFORM predic- 
tions for the numbers of Lya emitters at z > 7 were also given 
in Le Delliou et al (2006). We can see that the predictions from 
the different methods are fairly similar at z = 7.7, but grad- 
ually diverge from each other with increasing redshift. For the 
highest redshift, z = 12.5, the TM05 model fails in producing a 
prediction due to numerical problems. We note that making pre- 
dictions for z = 12.5 is challenging, for several reasons. Even 
though only ~ 200 Myrs separate the ages of the Universe be- 



4. 1. DaZle - Dark ages z Lyman-a Explorer 

DaZle is a visitor mode instrument placed on the VLT UT3 
(Horton et al. 2004). The instrument is designed to use a narrow- 
band differential imaging technique, i.e. observing the same field 
with two very narrow filters with slightly offset central wave- 
length. Objects with Lya in one of the filters can then be selected 
from the differential image of both filters. The field-of-view of 
DaZle is 6.83' x 6.83' and it is expected to reach a flux level of 
2 X 10~^* erg s~^ cm~^ {5a) in 10 hours of integration in one 
filter. This corresponds to a luminosity limit at redshift z = 7.7 
of log (LLyct) = 42.13 erg s^^. The two initial filters are centred 
on ZLya = 7.68 and 7.73 (with widths Az = 0.006 and 0.025 
respectively) and at this redshift, the surveyed volume becomes 
1340 Mpc'^ per pointing per filter pair. Thus, from Fig.|2] we can 
conclude that DaZle will discover ~ 0.16 — 0.45 candidates at 
z = 7.7 with one pointing and filter pair. 

4.2. ELVIS - Emission Line galaxies with VISTA Survey 

ELVISQ is part of Ultra- VISTA, a future ESO Public Survey with 
VISTAQ (Visible and Infrared Survey Telescope for Astronomy). 
Ultra- VISTA is planned to do very deep near-infrared broad- and 
narrow-band imaging in the COSMOS field. It will observe four 
strips with a total area of 0.9 deg^. The narrow-band filter is 
focused on the z^ya ~ 8.8 sky background window with cen- 
tral wavelength Ac = 1185 nm, and redshift width Az = 0.1. 
The flux limit of the narrow-band images is expected to reach 
3.7 X 10^^* erg s^^ cm^^ (5a) after the full survey has been 
completed. Ultra- VISTA will run from early 2008 for about 5 
years and all the data will be public. ELVIS is presented further 

' www.astro.ku.dk/^kim/ELVIS.html 
2 www.vista.ac.uk 
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Fig. 2. Predicted Lya luminosity functions at z > 7. Red lines are extrapolations from observed luminosity functions at lower 
redshift, green lines are TM05 models and black lines are GALFORM models. Different linestyles show different redshifts z = 7.7, 
8.8 and 12.5. No prediction is shown for the TM05 model at z = 12.5. Hatched area marks observational upper limits from Willis & 
Courbin (2005) and Cuby et al. (2006), both at redshift z = 8.8. In the left panel, the GALFORM predictions are shown for escape 
fraction /esc = 0.02 (our standard value), while in the right panel, they are shown for f^gc = 0.2. The predictions from the other 
two methods are identical in both panels. 



in Nilsson et al. (2006b). ELVIS will survey several different 
emission-lines (e.g. Ha at redshift z = 0.8, [OIII] at redshift 
z = 1.4 and [Oil] at redshift z = 2.2) as well as the Lya line. 

When the survey is complete, the final mosaic will reach a 
Lya luminosity limit of log (L^yct) = 42.53 erg s~^. The vol- 
ume surveyed will be 5.41 x 10'"' Mpc'^. From Fig.|2]we see that 
ELVIS should be expected to detect 3-20 LEGOs at z = 8.8. 

4.3. JWST 

A possibility even further into the future is to use the James 
Webb Space Telescop^l (JWST). JWST is scheduled for launch 
in 2013 and will have excellent capabilities within the near- and 
mid-infrared regions of the spectrum. Two of the instruments 
aboard JWST could be used for narrow-band surveys; NIRCam, 
the near-infrared camera, and TFI, the tunable filter imager (for 
a review on JWST see Gardner et al. 2006). NIRCam will have 
31 filters, of which nine are narrow -band filters. The filter with 
shortest wavelength has central wavelength Ac = 1.644 /im 
(F164N; ZLya = 12.5, Az = 0.135). TFI will have tunable fil- 
ters with variable central wavelength, however it is only sensitive 
at wavelengths larger than A ~ 1.6 /im. NIRCam is expected to 
reach a flux Umit of ^ Ix 10^^® erg s^^ cm^^ (5cr) in 10000 s of 
exposure time. Hence, a flux Umit of 5 x 10"^^ erg s~^ cm~^ 
(5(7, log (LLya{z = 12.5)) = 42.00 erg s~^) could be reached 
in 10 hours, assuming that the sensitivity is proportional to the 
square root of the exposure time. TFI is expected to be able 
to reach a flux limit almost a factor of two deeper in the same 
time, however it has a field-of-view of only half of the NIRCam 
(which is 2 X 2.16' x 2.16'). In one NIRCam pointing at redshift 
z = 12.5, approximately 1640 Mpc'^ are surveyed. Again, from 
Fig. 121 we can estimate that we will detect < 0.1 galaxies per 
10-hour pointing with NIRCam. However, the number of detec- 
tions depends strongly on the escape fraction which is unknown 
at such high redshifts, and thus the number of detected galaxies 
can be larger 
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5. Constraints on the early Universe 

Of the surveys at these redshifts that have been presented in pre- 
vious articles (Horton et al. 2004; Willis & Courbin 2005; Cuby 
et al. 2006; Nilsson et al. 2006b), or are conceivable (JWST, see 
section 1431 ) only ELVIS will detect a large enough sample to 
start to measure the luminosity functions and the extent of reion- 
isation at these redshifts and to study the fraction of PopIII stars 
in the population. We here discuss these issues with respect to 
ELVIS. 

From the semi-analytical modelling, we can make mock ob- 
servations of the ELVIS survey. The procedure to produce these 
catalogues is explained in detail in Orsi et al. (in prep.), but 
the outline of the process is that galaxies from GALFORM are 
placed in matching dark matter haloes in the Millenium N-body 
simulation (Springel et al. 2005), which is a cubical volume in 
a CDM universe of comoving size 500 Mpc/h, thus creating a 
mock universe with simulated galaxies which includes all the ef- 
fects of clustering. We can then make mock observations of this 
simulated Universe, including the same limits on flux, redshift, 
sky area etc. as for any real survey, and from these observations 
produce mock galaxy catalogues. From the mock catalogues, we 
can in turn make mock luminosity functions of Lya emitters at 
redshift z ~ 8.8. In Fig. [3] we plot the "observed" luminosity 
functions in the 112 mock catalogues taken from different re- 
gions of the Millenium simulation volume. Note that for mak- 
ing these mock catalogues, GALFORM was run with the same 
cosmological parameters as in the Millenium simulation itself, 
which are slightly different from the "concordance" values as- 
sumed elsewhere in this paper, as described in the Introduction 
(this is why the mean luminosity function for the whole sim- 
ulation volume which is plotted in Fig. |3] is slightly different 
from the GALFORM prediction for z = 8.8 plotted in Fig.©. 
We used escape fraction f^sc = 0.02. The figure shows that 
the spread in number density between the different mock cat- 
alogues is large, almost a factor of ten in number density in each 
luminosity bin. This is a consequence both of the small num- 
bers of galaxies in the mock surveys and of galaxy clustering. 
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log(LLyc,[erg s"']) 

Fig. 3. Luminosity functions at z = 8.8 for a set of mock ELVIS 
surveys computed using GALFORM. The 112 mock surveys 
are identical apart from being taken from different regions in 
the Millennium simulation volume. The open circles show num- 
ber counts in each mock catalogue, in four luminosity bins. The 
black dot with error bars shows the median of the mocks in each 
bin, with the error bars showing the 10-90% range. The thin lines 
are best fit Schechter functions to the median points with dif- 
ferent assumed faint end slopes. The thick solid line shows the 
"true" luminosity function, as measured from model galaxies in 
the total Millenium simulation volume. 



which causes "cosmic variance" between different sample vol- 
umes. The prediction from GALFORM is therefore that it will 
be difficult to accurately measure the luminosity function of Lya 
emitters at z = 8.8 even using the sample from the large area 
ELVIS survey. In particular, there will be no useful constraint 
on the faint-end slope a. This is simply a consequence of the 
flux limit of narrow-band surveys, i.e. even if we use the me- 
dian values of the luminosity function from all the mocks, then 
Schechter functions with slopes in the range —1 to —2 all give 
acceptable fits, as illustrated in Fig. [3] However, if all the data are 
combined in one luminosity bin, it should be possible to measure 
0* assuming values for a and L*. The possibility of including 
data points from several surveys at different luminosities (e.g. 
also lensing surveys that probe the faint end of the luminosity 
function) would also significantly improve the results. 

Two suggested methods of constraining reionisation from 
observations of Lya emitters, without requiring spectroscopy, 
are to measure the clustering of Lya-sources and to compare 
the Lya and UV continuum luminosity functions at these red- 
shifts (Kashikawa et al. 2006; Dijkstra, Wyithe & Haiman 2007; 
McQuinn et al. 2007). McQuinn et al. (2007) show that large 
Hll bubbles may exist during reionisation, and that these will 
enhance the observed clustering of Lya emitters in proportion 
to the fraction of neutral hydrogen in the Universe. A sample of 
~ 50 emitters will be enough to constrain the level of reionisa- 
tion using this effect (McQuinn, priv. communication), almost 
within reach of the ELVIS survey. A future, extended version of 



ELVIS would be able to place very tight constraints on reionisa- 
tion. In Kashikawa et al. (2006) and Dijkstra, Wyithe & Haiman 
(2007) the use of the combination of the UV and Lya LFs to 
constrain the IGM transmission is explored. Lya emission will 
be much more susceptible to IGM absorption than the continuum 
emission and thus the ratio between the two LFs will give infor- 
mation on the level of IGM ionisation. However, with increasing 
redshift for Lya, the continuum emission will be increasingly 
difficult to observe, and it is unclear if this method will be feasi- 
ble for surveys such as ELVIS. 

It is possible that galaxies at z = 8.8 still have a signifi- 
cant population of primordial PopIII stars. A test for the frac- 
tion of primordial stars is the amount of Hell 1640 A emis- 
sion (Schaerer 2003; Tumlinson, Schull & Venkatesan 2003). 
Depending on models, these authors predict that the Hell 1640 A 
emission line should have a flux between 1 — 10 % of the flux in 
the Lya line. For ELVIS z = 8.8 Lya emitters, the Hell 1640 A 
line is redshifted to 1.61 ^m. Due to the many OH sky emission 
lines in this region of the spectrum, it would be desirable to try 
to observe the Hell 1640 A line from a space-based observatory 
such as JWST. According to the JWST homepage, NIRSpec will 
achieve a sensitivity in the medium resolution mode on an emis- 
sion line at 1.6 fim of 7 x 10^^^ erg s^^ cm^^ (lOcr) for an 
exposure time of 10^ s (30 hours). Thus, if the Lya emission line 
has a flux of ~ 5 x 10~^^ erg s"^ cm"^, the Hell 1640 A will 
be marginally detected with JWST in 30 hours of integration, de- 
pending on the ratio of Hell 1640 A to Lya flux. The NIRSpec 
sensitivity increases at longer wavelengths, but the increasing lu- 
minosity distance to galaxy candidates with Hell 1640 A emis- 
sion at longer wavelengths will most likely counteract this effect. 

6. Discussion 

We summarise our predictions for number counts of Lya emit- 
ters in narrow-band surveys in Fig. H) We also summarise the 
numbers of detected objects for specific current and future sur- 
veys in Table [3] A few comments can be made on the differ- 
ences in predictions between the two models. Firstly, as can be 
seen in Fig. |4] and also Fig. [21 the luminosity functions have 
steeper faint-end slopes in the GALFORM models than in the 
TM05 models. Secondly, the GALFORM and TM05 models 
predict similar amounts of evolution at a given flux over the 
range z = 6 — 9 where they can be compared. 

Several factors enter into the error bars of our predictions. 
One problem is the uncertainties in, and disagreement between, 
the observed lower redshift luminosity functions which are used 
to calibrate the theoretical models. There are many caveats in 
producing Lya luminosity functions, of which the selection 
function is the most difficult to correct for The problem arises 
from that the filter transmission curve is not box-shaped, but 
rather gaussian. Thus, only brighter objects will be observed at 
the wings of the filter, and these will be observed to have smaller 
than intrinsic luminosities. Secondly, the equivalent width (EW) 
limit that the survey is complete to depends on the depth of the 
broad-band images used for the selection. Thirdly, if the sam- 
ple is a photometric sample, it is possible that there are lower 
redshift interlopers, where the emission line is e.g. [Oil], in the 
sample. Finally, the samples are still so small that we have to 
deal with small number statistics. All of these problems cause 
the observed luminosity function at lower redshifts to be uncer- 
tain. 

Both theoretical models (semi-analytical and phenomeno- 
logical) have uncertainties resulting from how they model the 
galaxy formation process, and also from the assumption that the 
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Table 3. Number of predicted/observed objects per observed field in several present and future surveys from two theoretical models. 
Data from Subaru XMM Deep Field (SXDS) are from Ouchi et al. (2005), Shimasaku et al. (2006) and Kashikawa et al. (2006). 
GALFORM predictions are made assuming an escape fraction of fesc = 0.02. 



Name 


Redsliift 


Area (arcmin'^) 


Luminosity limit {5a, erg s ^) 


GALFORM 


TM05 


Observed number 


SXDS (Ouclii) 


5.7 


8100 




443 


339 


515 


SXDS (Sliimasaku) 


5.7 


775 


■|^q42.40 


112 


86 


83 


SXDS (Kasliikawa) 


6.5 


918 


^q42.27 


108 


57 


58 


DaZle 


7.7 


47 


^q42.13 


0.45 


0.16 




ELVIS 


8.8 


3240 


■|^g42.50 


20 


2.8 




CubyOe 


8.8 


31 


^043.10 


0.0003 


0.0 





W&C05 


8.8 


6.3 


■|^q42.25 


0.015 


0.002 





JWST 


12.5 


9.3 


■|^q42.00 


0.018 
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Fig. 4. Summary of predictions. The plot shows the number of 
Lya emitting galaxies expected per square degree per redshift 
interval Az = 0.1 as a function of redshift and observed flux 
limit. The predictions of GALFORM are shown in black and of 
the TM05 model in red. The different line styles are for different 
flux limits. 



fraction of Lya photons escaping from galaxies is constant and 
does not change with redshift. In addition, neither model in- 
cludes attenuation of the Lya flux due to neutral hydrogen in 
the IGM. This attenuation would be expected to be strong at 
z > Zreion, whcu the IGM is neutral, and weaker at z < Zreion, 
when most of the IGM is ionised. The degree of attenuation de- 
pends on a number of different effects, as analysed in Santos 
(2004), and discussed in Le Delliou et al. (2006), and is currently 
very uncertain. Nonetheless, this attenuation is expected to pro- 
duce observable effects on the evolution of the Lya luminosity 
function, if reionisation occurs within the redshift range covered 
by future observations, and so estimating the reionisation red- 



shift and the neutral fraction after reionisation are included in 
the science goals of these surveys. 

It is apparent that the key to acquiring a large sample of Lya- 
emitting galaxies at redshifts greater than 7 is both depth and 
area. In a recent paper. Stark, Loeb & Ellis (2007) suggest that 
one of the most efficient means of finding very high redshift Lya 
emitters is through spectroscopic surveys focused on gravita- 
tional lensing clusters. Lensing surveys could easily reach down 
to a luminosity limit of lO**"'^ erg s^^ in a few tens of hours. 
However, the surveyed volumes are very small, of the order of a 
hundred Mpc^^. For a lensed survey, the area in the source plane 
is reduced by the same factor that the flux is amplified, so in 
principle one gains in the total number of objects detected rel- 
ative to an unlensed survey if the luminosity function is steeper 
than N{> L) oc L-'^. In the GALFORM and TM05 models, the 
asymptotic faint-end slope is shallower than this, but at higher 
luminosities, the slope can be steeper. For example, GALFORM 
predicts that at z = 10, the average slope in the luminosity range 
IQ-H-io-is erg s^i is close to iV(> L) cx L"^ (ggg pjg^ g 
Le Delliou et al. 2006), so that a lensing amplification of 10 re- 
sults in 10 times more objects being detected, with intrinsic lu- 
minosities 10 times lower, compared to an unlensed survey with 
the same area and flux limit. Therefore lensing and narrow-band 
surveys are complementary to each other as they probe differ- 
ent parts of the luminosity function. With either type of survey, 
reaching a significant sample of redshift z ^ 7 — 8 should be 
possible in the next few years with telescopes/instruments in use 
or soon available. 

An interesting type of object found recently in narrow-band 
surveys are the Lya blobs, large nebulae with diameters up to 
150 kpc and Lya luminosities up to 10^^ erg s^^ with or with- 
out counterpart galaxies (e.g. Steidel et al. 2000; Matsuda et 
al. 2004; Nilsson et al. 2006a). Several mechanisms have been 
proposed to explain this phenomenon, including starburst galax- 
ies and superwinds, AGN activity or cold accretion. It is inter- 
esting to consider if such objects would be detected in any of 
these surveys, assuming they exist at these redshifts. A typical 
Lya blob will have a luminosity of ^ 10"^^ erg s^^ and a ra- 
dius of, say, 25 kpc. This will result in a surface brightness of 
~ 5 X IQ-^^ erg s~^ kpc^'^. Thus, a narrow-band survey will 
have to reach a flux limit, as measured in a 2" radius aperture 
of 1.3 X 10**^ erg s^^ at redshift z = 8.8, corresponding to 
logL = 42.11. (An aperture radius of around 2" is expected to 
be roughly optimal for signal-to-noise.) For lower or higher red- 
shifts, this limit is higher or lower respectively. Thus, ELVIS will 
not be able to detect Lya blobs unless they are brighter and/or 
more compact at higher redshift than a typical blob at lower red- 
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shift. DaZle and JWST could in principle detect this type of ob- 
ject, but only if they are very abundant in the very high redshift 
Universe, due to the small survey volumes of these instruments. 
It is of course highly uncertain what properties such Lya blobs 
would have at z ~ 7 — 9, or their space density, but it appears 
unlikely that the future surveys presented here would detect any 
such objects. 

To find compact Lya emitters at redshifts z > 10 in signif- 
icant numbers we will probably have to await instruments even 
further in the future. If a future 40-m ELT (Extremely Large 
Telescope) was equipped with a wide-field NIR imager and a 
narrow-band filter of similar width to ELVIS, it could reach a 
luminosity limit of L lO**^'^ erg s~^ at redshift z = 10.1 
(where a suitably large atmospheric window exists) in approxi- 
mately 20 hours. Using the GALFORM model for z = 10, the 
number density should be N(>L)»i 4 x 10^"^ Mpc^'^ at this lu- 
minosity limit. Thus, to get a sample of ten Lya emitters would 
require imaging an area on the sky of approximately 16 square 
arcminutes, assuming a narrow-band filter with redshift range 
10.05 < ZLya < 10.15. This could be achieved with one point- 
ing if the detector has a field-of-view of 6 arcmin on a side, as 
suggested by the ESO ELT Working GroupQ. It should of course 
be noted that these are very tentative numbers, but they display 
the possibilities of far future instruments. 
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